Optimal tracking for pairs of qubit states 
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In classical control theory, tracking refers to the ability to perform measurements and feedback on 
a classical system in order to enforce some desired dynamics. In this paper we investigate a simple 
version of quantum tracking, namely, we look at how to optimally transform the state of a single 
qubit into a given target state, when the system can be prepared in two different ways, and the 
target state depends on the choice of preparation. We propose a tracking strategy that is proved to 
be optimal for any input and target states. Applications in the context of state discrimination, state 
purification, state stabilization and state-dependent quantum cloning are presented, where existing 
optimality results are recovered and extended. 

PACS numbers: 03.67.Pp, 03.65.Ta, 03.67.-a 



I. INTRODUCTION 



A common goal of many problems in quantum infor- 
mation science is the search for quantum operations that 
simultaneously transform a set of given input quantum 
states into another pre-specified set. Well known exam- 
ples are task such as quantum cloning, state discrimina- 
tion and quantum error correction. 

In general, though, quantum mechanics forbids arbi- 
trary quantum state dynamics. As a result, one is left 
with several examples of "impossible quantum machines" 
0. Not only is quantum cloning unachievable j^, ||, 
but also quantum state discrimination strategies are typ- 
ically subject to non-zero misidentification probabilities 
H and/or inconclusive outcomes 0, ||] and there are 
no quantum error correction protocols capable of fully 
reverting the action of an arbitrary noise model B . 

Nevertheless, it is still possible to approximate ideal 
(but unphysical) transformations with optimal (but 
physical) ones. This provides quantum limits to the per- 
formance of tasks such as state discrimination, cloning 
and so on. In this paper, we study the general problem 
of transforming the state of a single qubit into a given 
target state, when the system can be prepared in two dif- 
ferent ways, and the target state depends on the choice 
of preparation. We call this task quantum tracking, a 
term borrowed from classical control theory. Our main 
result is an analytical description of an optimal quantum 
tracking strategy. 

More specifically, the quantum tracking problem stud- 
ied here can be understood as follows. Consider that 
Alice prepares either a qubit state p\ with probability 
7Ti or p2 with probability TT2- Bob is allowed to interact 
with the system in any physically allowed way, aiming to 



enforce the tracking rule 



if Alice prepared pi , then ouput p i 



(1) 
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for i = l,2 and some given qubit density matrices ~p i . 

At his disposal, Bob has all the information about the 
possible preparations pi and their respective prior prob- 
abilities 7Tj, but not the actual preparation (the value of 
the index i). 

Because quantum states are generally not perfectly dis- 
tinguishable, a strategy that attempts to identify Alice's 
preparation and then reprepare the target according to 
rule (Q) is not always guaranteed to succeed. In fact, this 
limited distinguishability is an unsurpassable obstacle in 
the implementation of ([!]). 

In this paper, an optimal solution is defined as follows. 
Amongst all the physical transformations acting on the 
input states pi, an optimal one is any map that out- 
puts density matrices p\ such that the averaged Hilbert- 
Schmidt inner product between p[ and 'p i is maximal. 
When 'p i are pure states, such a figure-of-merit coincides 
with the averaged Uhlmann-Jozsa fidelity jnj |Tl]| , and 
this notion of optimality gains an appealing operational 
interpretation ]l2] ]. Suppose that Alice [aware of her 
preparation and of rule ([[])], decides to check whether 
Bob prepared the density matrix he was supposed to, and 
for that purpose she performs a verification measurement 
on the density matrix produced by him. If Bob chooses 
an optimal transformation according to the above pre- 
scription, then the probability he will pass Alice's test is 
as large as allowed by quantum mechanics. 

The tracking problem resembles the transformability 
problem for pairs of qubit states studied by Alberti and 
Uhlmann in the 80's (see also Appendix |a|). In 
a criterion based on the distinguishability between the 
source density matrices and the distinguishability be- 
tween the target density matrices was developed in order 
to decide on the existence of a completely positive and 
trace preserving (CPTP) map simultaneously transform- 
ing each source into each target. 

Although Alberti and Uhlmann's criterion classifies the 
set of states pi, p2, ~p~\ and ~p 2 f° r which rule (|j) can be 
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The paper is structured as follows. Section [Ej] intro- 
duces the formal statement of the problem and our work- 
ing strateg y, w hich is proved to be optimal in Section III 
In section [V we evaluate the performance of the opti- 



mal strategy in the contexts of quantum state discrimi- 
nation, quantum state stabilization in the presence noise, 
perfect quantum tracking, and state-dependent quantum 
cloning. Section [v] proposes a physical implementation 
of our strategy in terms of closed and open loop control. 
Section [vj discusses generalizations of the problem and 
concludes the paper. 



II. PROBLEM AND STRATEGY 



FIG. 1: (Color online) The transformation of the mixed states 
p\ and p2 (in blue) into the pure states ~p 1 and p 2 (in red) is not 
a physical one. However, there is a physical transformation C 
capable of transforming the input states into states C(pi) (in 
black) which closely approximates the targets. Note, in the 
detail, that C(pi) and C(p2) are still slightly mixed. 



satisfied, it does not provide a construction of the CPTP 
map implementing that transformation, nor touch the 
problem of how to find a feasible approximation when 
the criterion is not satisfied. For many purposes, the re- 
quirement of perfectly converting sources into targets is 
unnecessarily strong, as some strictly impossible physical 
transformations can be very well approximated by physi- 
cal ones, as illustrated in Fig[l]. In fact, any experimental 
realization of a map is just an approximation of it. 

Another problem closely related to the aims of this 
paper was investigated in Ref. |l4|j . Specifically, we con- 
sidered the problem of determining the optimal quantum 
operation to stabilize the state of a single qubit, ran- 
domly prepared in one of two pure states, against the 
effect of dephasing noise. The results of are here ex- 
tended in several ways. The input states are allowed to 
be mixed and prepared with arbitrary prior probability 
distribution; the noise model is arbitrary and, most im- 
portantly, the stabilization task is replaced with tracking. 

Finally, there is an intrinsic connection between the 
quantum tracking problem and the "optimization ap- 
proach" PLM HI' HI' Hi t0 q uantum error 
correction ]41|. In these references, the encoding and re- 
covery operations are regarded as optimization variables 
whose optimal values maximize a given figure-of-merit 
(typically a function bounded between and 1, equal to 
1 if and only if the noise dynamics is reversible). Efficient 
numerical methods are then proposed to solve the opti- 
mization problem. The key differences between our work 
and these references is that we do not consider encoding 
of the initial state and focus on reverting the noise dy- 
namics experienced by only a pair of states. By doing so, 
the optimization of the recovery operation can be handled 
analytically for a conveniently chosen figure-of-merit. 



In this section we give a formal statement of the prob- 
lem of interest and introduce our strategy 

A. The Problem 

Formally, the problem we set out to solve can be stated 

as 

Problem. Given qubit density matrices p\, p 2 , ~p 1 , ~p 2 
(with pi ^ P2) o,nd probabilities 7Ti, tt 2 with 7Ti + 7T2 = 1, 
find a quantum operation C maximizing 

^hs = TTiTr [C(/>i)pJ + ^ 2 Tr [C(p 2 )p 2 ] . (2) 

We will refer to this as "the tracking problem" . 

The choice of the average Hilbert-Schmidt inner prod- 
uct as our figure-of-merit J-rs is motivated by technical 
reasons (to be clarified later), and by the fact that for 
pure target states (the case of greater interest as far as 
applications are concerned), ^hs is precisely equal to the 
average fidelity. When the target states are mixed, .Fhs 
is a lower bound to the average fidelity Although 
not as well motivated as in the case of pure target states, 
the determination of the quantum operation C maximiz- 
ing JFhs can still be useful for mixed target states. For 
example, if a certain application requires tracking to be 
performed with average fidelity / and the optimal value 
of our figure-of-merit is such that J-"hs > f , then C is 
suitable for the task. 

As a final remark, note that we do not exclude the 
case ~p 1 = ~p~2 from the statement of the problem. How- 
ever, we will exclude the case 'p 1 = ~p 2 = 1/2 from the 
following analysis. Obviously, this particular transforma- 
tion is always feasible and achieved with the completely 
depolarizing channel. 

Next, we propose a strategy that will later be proved 
to be a solution of this tracking problem. 

B. The Strategy 

In this section, we provide an analytical solution of the 
tracking problem, i.e., we detail the structure of an opti- 
mal tracking operation C and derive closed forms for the 
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associated maximal value of the figure-of- merit J-hs- The 
scheme proposed here was constructed by incorporating 
some features observed from the numerical solution of 
the tracking problem in an analytical optimization pro- 
cedure. In the next section, we will show that the track- 
ing problem can be cast as a semidefinite program (SDP) 
p2| , [2^ |, and will employ the theory for this type of op- 
timization problem to prove that the strategy presented 
here actually solves the tracking problem. 

A quantum operation is a description of a certain phys- 
ically allowed evolution of a quantum state. For a closed 
quantum system (not interacting with an environment) 
this description is given by the familiar unitary evolution 
of Schrddinger's equation. For open quantum systems, 
unitary evolution alone does not account for every possi- 
ble state transformation — in this case, the set of quan- 
tum operations is identified with the more comprehensive 
set of completely positive and trace preserving (CPTP) 
maps. 

Any one qubit CPTP map C can be decomposed as 
@@ 

C(?) - UV(V^)U^ , (3) 

where U, V are unitary matrices and T> induces an afiinc 
transformation on the input Bloch vectors; namely, it 
contracts the x, y and z components via a multiplicative 
factor, and subsequently adds a fixed number to them. 
Any non-unitary evolution arises from a transformation 
of this type. In the framework of Eq. (^), unitary dy- 
namics is simply obtained by making the affine map T> 
redundant (e.g., multiplying by l's and adding O's to the 
x— , y— and z— Bloch components). 

In general, CPTP maps reduce the distinguishability 
of quantum states. On the Bloch sphere this typically 
corresponds to a reduction of the Bloch vector length 
and angles between vectors. In contrast unitary dynam- 
ics preserves the angles between Bloch vectors and their 
lengths. For the tracking problem, we can imagine that 
in some cases the optimal strategy will preserve lengths 
and angles, i.e. it will be some unitary correction. We 
will construct an "indicator function" which will flag this 
case. 



1. Indicator function 

To gain some intuition, we start by constructing an 
indicator function for the simplest case of tracking with 
uniform priorities tt\ = ~ki — 1/2 from pure states (pi 
and P2) to pure states (p 1 and p 2 )- Throughout, and 
O will denote the angles between the Bloch vectors of 
pi, P2 and P27 respectively. It will also be convenient 
to define 9 and to be the half-angle between the Bloch 
vectors, i.e., 29 = and 28 — 0. 

Given that all the states involved are pure, it is 
straightforward to conclude that if = 0, then unitary 
dynamics is the best choice — a suitable rotation of the 
Bloch vectors of the inputs can perfectly bring them to 



coincide with the Bloch vectors of the targets (as opposed 
to a non-unitary evolution that would decrease the angle, 
hence excluding the possibility of perfect tracking) . 

A corollary of a theorem by Alberti and Uhlmann |ll| 
(see appendix ^) implies that any pure state transforma- 
tion such that > 0, can be perfectly implemented. If 
that is the case, then this transformation must be non- 
unitary, since a unitary would not be able to bring the 
angles to perfectly match. This suggests the introduction 
of the function 

Q:=0-0, (4) 

to indicate non-unitary dynamics whenever 51 > 0. Next, 
we argue that Q < indicates unitary dynamics, thus 
establishing ft as an example of indicator function we 
were looking for. 

We have already seen that Cl = implies unitary dy- 
namics. Intuitively, this conclusion can be extended to 
£1 < with the following reasoning. If < 0, any fur- 
ther decrease of the initial angle can only further separate 
the resulting states from the targets. Since there is not 
any quantum operation capable of increasing this angle, 
the best policy must be to preserve it, hence a unitary. 

The above discussion may suggest that the optimal 
indicator function is merely a comparison of the distin- 
guishabilities between sources and targets. If the sources 
are more distinguishable than the targets, then we em- 
ploy a quantum measurement to decrease the distin- 
guishability, hence approximating the targets. If the 
sources are no more distinguishable than the targets, 
then we employ a unitary operation to avoid a further 
decrease of the overlap between the output states and 
the targets. Although this reasoning is certainly in agree- 
ment with the indicator function introduced above for the 
special case of pure states, it does not extend to mixed 
state transformations fi^ . 

If the states are not pure and the priorities are not 
uniform, it is much more difficult to understand how pu- 
rities, angles and priorities combine to form a meaningful 
decision criterion about the nature of the best dynamics. 
In order to introduce an indicator function for this gen- 
eral case (obtained from some mathematical optimization 
procedure, not from an heuristic argument), we first de- 
fine some useful notation. 

Let Ri and Ri be the Bloch vectors for pi and ^{pi , 
respectively (note that Ri is not normalized). Symboli- 
cally, for 11 e {R,R}, define 

K+ := Tlx + TZ 2 , (5a) 
TZ := ni-112, (5b) 
TZ X := 72.1 X K2 , (5c) 

and as usual, the corresponding unbolded type gives the 
Euclidean norm 7?.{ + ,_. x } = ||7?.{ +i _ iX }||. Also, the fol- 
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lowing will be important throughout 

2 

T:= (l-Ri'Rj)(Ri-Rj), (6) 



S:= yJT 2 +4R X {R 2 _- Rl 



(7) 



We note tha t R i ^ i?2 guar antees that R 2 _ — R 2 X > 
(see Appendix B 1, Lemma B.l), hence both S and T are 
real numbers. 

In terms of these quantities, we define 



il:=S + T-2R x R > 



(8) 



with Q > indicating that non-unitary dynamics are 
required (which will be detailed as "procedure A") and 
f2 < indicating that unitary dynamics ( "procedure B" ) 
are required. 

Although it would be difficult to motivate the indicator 
function fl of Eq. (||) as we did with in Eq. (Q), it 
is possible to see that the former is equivalent to the 
latter in the case of pure qubit states. This is shown 
in Fig. ^, where it is also noticeable that even a simple 
generalization of the input states from pure to mixed 
states with the same level of mixedness (as measured by 
the norm of their Bloch vector R), is already sufficient 
to give a fairly non-trivial division line between the two 
types of dynamics. 



2. Procedure A 

In this section we present the details of the map C 
from Eq. (^) for > (which indicates non-unitary 
dynamics). 

Step 1. The rotation by the unitary V takes the two 
input Bloch vectors to vectors R[ and R' 2 in the xz- 
plane in such a way that they share a common positive 
a;-component and R[ ■ z > R' 2 ■ z, explicitly 



Ri 



TV 

It, = X 

1 R- 



(Ri-R-) 
R 



(9) 



Step 2. The afhne transformation T> shortens the x, 
y and z components of its inputs by multiplying them, 
respectively, by fit, ^2, M3 with < ^{1,2.3} < 1 an d 
subsequently adding s\ to the x component. Applied to 
RL that reads 



R' 



R'! 



(Ri-R-) , , 
V l R > z. (10) 



That such a transformation can be physically imple- 
mented is not a trivial fact. Indeed, strict conditions 
involving the parameters /i{i,2. 3} and s\ must be satis- 
fied to guarantee the feasibility of transformation (|ITj) as 
a CPTP map |^5l . The following values can be shown to 
satisfy these conditions 



0.75 




0.25 



0.25 0.5 0.75 

Source fidelity 



FIG. 2: Lines separating unitary and non-unitary dynamics, 
as prescribed by the indicator function of Eq. . For source 
states, we consider pairs of mixed states with Bloch vector 
length R and separated by an angle Q — 28 (Bloch sphere 
angle). For target states, we consider pairs of pure states 
separated by an angle = 29. The fidelity between the 
source states (horizontal axis) is 1 — R 2 sin 2 6, and the fidelity 
between the target states (vertical axis) is cos 2 8. The region 
where unitary dynamics is advisable (Q < 0) is indicated 
with an arrow. The intuitive notion that a measurement is 
employed when the targets are less distinguishable than the 
sources (and a unitary, otherwise), only holds if R — 1 (pure 
sources). For R = 0.9,0.8 the division line moves down in 
such a way as to increase the portion of the parameter space 
where non-unitary dynamics is advisable (fi > 0). 



Mi = 2^ 



;R„RxR- 



S(S + T) 3 x x ' 
2 



/'2 



M3 



Si = 



S + T 



RxRx , 



R X R- , 



S{S + T) 



2S(S + T) S 



{S + Tf - AR y R\ 



(11a) 
(lib) 
(11c) 
(lid) 



In Appendix [B l| we show that the only circumstances 
under which the inequalities S > and S + T > are 
not simultaneously satisfied have Q = 0. Therefore, the 
quantities above are real and well-defined for the present 
procedure (SI > 0). It is not difficult to check that /ii = 

M2M3 arL< i s i = — Mi)(l — M3)? so the resulting map 
acting on density matrices is an extremal point of the 
convex set of CPTP maps [p5l . 
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Remarkably, if the target Bloch vectors R\ and R 2 are 
parallel or anti-parallel (i.e., R x = 0), Eqs. (|TT| ) simplify 
to fii = fi 2 = M3 = and s\ = 1. This implies that 
R'l = x for i = 1, 2, or equivalently, that I? outputs 
1+) = (1°) + |1))_/V2 independently of the input. 

S^ep 5. For R x ^ 0, the unitary U rotates the input 
vectors R" to lie on the plane determined by the target 
vectors, and within that plane by a suitable angle. For 
R x = 0, U simply rotates the vector x in order to align 
it with _R + . In either case, U can be expressed as the 
following map 



R'l R'l' = k a R t + k l2 R 2 , 



(12) 



with 
k i 



a 2 + ft&JJx + (-l)^a (/3i - 02) Rr ■ R 1 



(13) 



—2 



|/3 1 i2 1 +^ 2 -R 2 || 2 + 2a (ySi — ^2)] J2 



(14) 

where we have defined the involutory map ~ such that 
1 = 2 and 2 = 1; and 



a = R 'l ■ x 



IS + T 
2S 



Rx ys(s + T) 



Ri ■ R 



(15) 
(16) 



With this, the figure-of-merit of Eq. (|2|) can be shown 
to be 



T A - - 4- -i 



where T a is obtained by substituting Eqs. ( |l5| ) and Jl6| ) 
into Eq. (Hj) and reads 



S + T 



(18) 



Fig. |[hlustrates this sequence of transformations for 
the case R x ^ 0. 



3. Procedure B 

As pointed out before, if O < then the affine transfor- 
mation T> is not implemented and the product VU gives 
the unitary dynamics. In this case, V can be chosen pre- 
cisely as in step 1 of procedure A. 

For R x ^ 0, we preserve the form of the transforma- 
tion U from Eq. (|l2|), but the values of a and f3i are 
given by 

Q' = Rj ■ x = -77— , (19) 



13, 



R- ' 
R'i ■ z Ri ■ R^ 



R> 



R X R- 



(20) 



For i? x = 0, assume that the Bloch vectors R\ and R 2 
are anti-parallel (this is without loss of generality, since 
parallel targets always exhibit Q > p3[ ). In particu- 
lar, take R\ parallel to z and R 2 parallel to —z. The 
following transformation specifies U in this case 



ryl , u , r>n 
±1.4 I ' -tl j 



where 



R> 
R- 



cost) 



R^ 

R- 



sin 



(Ri ■ R ) 
R- 
(Ri ■ R- 
R~ 



x 



sin-d — 



R X (R 1 -R 2 ) 



(21) 



(22) 



and cost? = — sin 2 d. We note that 1) is a valid 

angle since the rhs of Eq. (^2|) is bounded between —1 
and 1_Q. 

If /?i and i?2 do not align along the z direction as 
specified above, we simply apply a further rotation that 
aligns the z axis with the direction Ri / R\ . 

In both R x 7^ and R x = cases, the average fidelity 
can be computed to be 



fB _ ■ 



(23) 



2 2 

where is obtained by substituting Eqs. ([[9]) and fl2C| ) 
into Eq. dbf). After some manipulation we find 



R\ 



T + 2R X R > 



(24) 



III. OPTIMALITY PROOF 



In this section, we employ duality theory for SDPs to 
prove the following theorem 

Theorem III.l. Our tracking strategy (as described in 
Sec. II E ), implements optimal tracking between any pair 



of source and target qubit states and is, 
lution of the tracking problem introduced in Sec. II A 



therefore, a so- 



In the subsequent proof of this theorem, some famil- 
iarity with SDP theory is assumed. Standard reviews on 
the topic are |?2| |23| . More closely related to our pur- 
poses is [ p6| , where the connection between optimization 
of quantum operations and SDPs was first noted. Also 
relevant is Ref. fll4| , where a similar technique was used 
to approach a particular case of tracking. 



A. The tracking problem as a SDP 

We start by showing that the tracking problem can be 
formulated as a SDP. Formally, it can be written as 



max Tr [C(pi)-Kipl\ . 

C£{CPTP}^ 



(25) 
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FIG. 3: (Color online) Bloch sphere schematics of procedure A. a) The source (target) Bloch vectors Ri and R2 (Ri and R2) 
determine the plane ir s (nt)- b) V implements the rotation transforming 7r s to the xz-plane, in such a way that the vector 
R'x — R'2 is parallel to +z. c) The map D deforms the Bloch sphere into an ellipsoid of semi-axis fii, p2 and /i3 and translates 
it by si along the x axis. The resulting ellipsoid touches the original Bloch sphere — a feature related to the fact that T) is 
an extremal CPTP map. d) U rotates the resulting states to the plane 7r t , and within that plane by some angle such that the 
resulting states R"' and R'2' approximate Ri and R2, respectively. 



It will be convenient to rewrite C{p t ) as [^7j 

C(pi) = Tr A [(pj ® 1)K C ] , 
where Kc is the (unnormalized) Choi matrix 
K C = (X®C)(|*+)(*+|), 



(26) 



(27) 



and |*+) = 1 00) + |11). Eqs. (|2J) and (g7|) establish 
a one-to-one relation between the set of CPTP maps on 
qubits and the set of (unnormalized) 2-qubit density ma- 
trices satisfying Tt b K c = 1 @, J|| |jj |(J . Here, Tr A{B) 
denotes the partial trace operation over the first (second) 
qubit. 

Using this isomorphism, a straightforward manipula- 
tion gives for the objective function in ( p5|) the form 
-Tr[F K c ], where 



i=l 



(28) 



whereas the constraint C € {CPTP} becomes Kc > 
and Trg [Kc] = 1a- In conclusion, the tracking problem 
assumes the standard form of a SDP 



maximize — Tr [i^ifg] 
subject to Kq > 

Tr [(tjfc ® t)K c ] 



28, 



k0 ■ 



-0,1,2,3. 

(29) 

A special feature of (|2^) which will be explored next is 
that the replacement p. L 1— > ApiA' and 'p i 1— > B^j^B, with 
A,Be U(2), yields another SDP (with F n replacing F ) 
that achieves exactly the same optimal value. This can 

be easily seen by noting that: (i) Tr [Fo-RTc] = Tr FqKc , 

for K c = (A T ® B)U< C (A T ® B) and (ii) K c satisfies the 
constraints in pfjl) if and only if Kc does. 



B. The duality trick 

In terms of the Choi matrix, the strategy described in 
[IB| is written as K c = (V T ® U)K V (V T ® U)\ 



section 
with 



- 1 ® 1 + — 1 <g> X + y X ® X - 



M2 



/';! 



(30) 

Given the reasoning of the previous section, our strategy 
constitutes an optimal solution to the tracking problem 
if and only if the following SDP is solved with K = Kx>, 



maximize — Tr FqK 
subject to K > 



Tr 



(o*®l)Jf =2tf fc0 , fc = 0,l,2,3. 



with 



(31) 
(32) 



The above SDP has the strong duality property, i.e., its 
optimal value is guaranteed to be identical to the op- 
timal value of its dual problem |23|]. This fact follows, 
for example, from the "strict feasibility" of the point 
Kf — 1 ® 1/2, which satisfies the constraints of ( |3~l| ) 

with the strict inequality Kf > 0. 

From duality theory for SDPs, the problem above is 
solved with K = Kt> if and only if (i) K-p satisfies the 
constraints of ( |3l| ) and (ii) the linear matrix inequality 

F = F + ? oa®a+yiA"®i+y 2 ^®a+y3^®i > o (33) 

is satisfied by some quadruple (po>?i>?2,y3) such that 



2?o 



Tr F K V 



(34) 



7 



If that is the case, then the so-called "complementary 
slackness" condition Kt>F = 0, holds for the appro- 
priate values of coefficients r , pi, y 2 and £3. 

To see that (i) is verified, recall that the values of /Ltj.,2,3 
and S\ were chosen to make of T> an (extreme) CPTP 
map. As mentioned before, the Choi matrix of any such 
map (on qubits) is characterized by the constraints of 
problem (pl|). 

is merely .F^s or 



F K V 



For (ii), first note that — Tr 

■^HS gi ven m Eqs. ( p"7j ) and (p3|), depending on whether 
f2 > or fi < 0. In our particular problem, the comple- 
mentary slackness condition results in sufficient indepen- 
dent linear equations that joj Fi , ?2 and £3 are defined 
precisely. We find y 2 — and 



• If SI > 0, 

Fo = ^(i + r Q ) , 
fl = ^L(i + r B ), 

r ; ■■ ■■ -777- (niRi + n 2 R 2 ) ■ R 



4i?_ 



(35a) 
(35b) 

(35c) 



IV. EXAMPLES 

In this section we evaluate our tracking strategy at 
work in some physically relevant problems such as quan- 
tum state discrimination, quantum state purification, 
stabilization of quantum states in the presence of noise 
and state-dependent quantum cloning. Moreover, we also 
discuss the application of our strategy in circumstances 
where tracking is known to be perfectly achievable. The 
analysis presented in this section is meant to give an ex- 
plicit account on the wide range of physical applications 
of the tracking problem and its optimal solution. 



A. Quantum State Discrimination 

A standard result in quantum state discrimination is 
the Helstrom measurement |j| , which consists of a projec- 
tive quantum measurement that maximizes the probabil- 
ity (Pfjoist) of correctly identifying the state of a quantum 
system that could have been prepared in two different 
states. Describing the possible preparations by p\ with 
probability p\ and p 2 with probability p 2 , the Helstrom 
measurement gives 



• If fi < 0, 

fo = J(i + r 6 ) , 



Heist 



PlPl -P2P2\W, 



(40) 



F3 



1 



4i?_ 



(tti.Ri + tt 2 R 2 ) ■ R + — 

1 b 



where, for brevity, we have defined 



S := £ (J* • R-) (Ri ■ R+) , 

1=1 

£ : = RxR+ + RxR 2 - , 
which can be shown to satisfy the relation 



S 2 +£ 2 = R 2 _R 2 + T 2 b . 



(36a) 
(36b) 

(36c) 



(37) 
(38) 

(39) 



We prove in Appendix B 2 that, although T a and T b ap- 
pear in the denominator of some of the coefficients in Eqs. 
( |35| ) and (|35|) , no singularities occur if the indicated range 
of fl is observed. 

With the set of coefficients (f|) and @, Eq. @ is 
clearly satisfied. As a result, the optimality of our track- 
ing strategy is solely dependent on proving the linear 
matrix inequ ality F > for the above set of coefficients. 
In Appendix B 3 we study the characteristic polynomial 



of F and conclude that all of its roots are non-negative, 



where || • ||t r denotes the trace norm. 

In this section, we propose a quantum state discrim- 
ination protocol for a pair of qubit states based on the 
tracking strategy introduced in the last section. We will 
show that it is equivalent to Helstrom's strategy, as it 
will give the same correct identification probability of 
Eq. ©. 

Our quantum state discrimination protocol consists of 
two simple steps: first we apply the tracking operation 
C introduced in section [II B| to approximate the states 
to be discriminated to some pair of orthogonal states. 
Without loss of generality, we take ~p x = |0)(0| and ~p 2 = 
|1)(1|. The priority of each transformation is taken to be 
identical to the prior probabilities with which p x and p 2 
are prepared, i.e., 7ri = Pi. As the second and final step, 
we perform the quantum measurement {|0)(0| , |1)(1|}, 
under the understanding that an outcome '0' suggests 
the preparation to be pi and an outcome T' suggests p 2 . 

The probability of a correct identification under this 
tracking scheme is given by Born's rule, averaged with 
the prior probabilities, 



Ptrack = PiTr [C(pi) |0)(0|] + P2 Ty [C(p 2 ) |1)(1| 
= 7nTr[C(pi)p 1 ]+^ 2 Tr[C(p 2 )p 2 ] . 



(41) 



thus proving Theorem III. 1 



By comparing Eqs. ( f4l| ) and (|^), one promptly recognizes 
that Ptrack is precisely the performance of the operation 
C for tracking from pi to pj with priority 7Tj, as measured 
by Fhs- Hence, in the case of pi = |0)(0|, p 2 = |1)(1| 
and TTi = pi, Eqs. (0) and (|2^) give the probability 
of success of our discrimination scheme for Q > and 
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51 < 0, respectively. Next, we make these formulas more 
explicit. 

Using the condition R x = in Eqs. (Q) and (||), we ob- 
tain fl = \T\ + T. Essentially, this means that T assumes 
the role of the indicator function: if T > 0, then > 
and we employ procedure A; if T < 0, then Q = and we 
employ procedure B. Substituting Ri = —(—lYpiZ into 
Eq. some simple algebra gives 



T=(pi 



and we can write 



P2) 



llpi-Ri -P2R2 



track — 



2 1 2 
\PlRl 



\\Pl -P2| 
P2-R2I! 



if T > 
if T < 



(42) 



(43) 



where the first line follows from Eq. fll7| ) and the second 
fromEq. @. 

It is a tedious exercise (essentially the computation of 
the eigenvalues of pipi — P2P2) to re-express Eq. (^0|) in 
terms of the Bloch vectors R4. The result is exactly 



Ft 



Heist 



Pi 



track ■ 



(44) 



hence establishing the claimed equivalence between our 
strategy and Helstrom's. 

Note that if T > 0, Ptrack is independent of the states 
we are trying to distinguish, but merely dependent on 
the probabilities with which they occur. This can be un- 
derstood by looking at the details of the affine operation 
taking place in procedure A. As noted before, for R x =0 
(as is the case for orthogonal targets), the affine map is 
such that /ii = fj,2 = fJ-3 = and s% — 1; that is, the 
source states are completely depolarized and a new state 
|+) is prepared instead. Next, this state is rotated by the 
unitary U and the measurement is finally performed. 

It is easy to see that for pi = p 2 = 1/2 the condition 
T > (procedure A) never holds. However, as we deviate 
from the uniform distribution, the volume of the param- 
eter space where procedure A is recommended grows to 
fully cover the space when pi = or p\ = 1. This is 
shown in Fig. . 



B. Quantum state Purification 

In this section, we consider a kind of state purification 
task where we aim to transform a pair of mixed source 
states pi and P2 with the same degree of mixedness (R\ = 
R2 = R < 1) that are separated in the Bloch sphere by an 
angle 29 € (0, tt] into a pair of pure target states ~p 1 and p 2 
separated by the same Bloch sphere angle 26. In other 
words, our purification task consists of elongating the 
Bloch vectors while preserving the angle between them. 

For later use, it will be convenient to derive formulas 
for the indicator function and figure-of-merit of a slightly 
more general problem, where the angle between the tar- 
get Bloch vectors is 28 £ [0,ir]. The purification task can 
be recovered by restricting 9 — 9. In addition, we will 




FIG. 4: (Color online) Each color of sheet represents a fixed 
deviation from the uniform probability distribution. The 
sheets divide the parameter space in two regions. The ar- 
rows designate the regions where procedure A (non-unitary) 
is recommended. On the other side, procedure B (unitary) 
is recommended. The larger the deviation from pi — 1/2, 
the larger the region where a non-unitary preparation for the 
measurement {|0}(0| , is advisable. 



allow the priorities of the transformations pi — > J> i to be 
arbitrary positive scalars -Ki such that 7Ti +7T2 = 1. Later, 
we make w± = 7T2, in order to simplify the formulas. 

In this generalized purification framework (9^9), the 
indicator function is obtained from Eq. (g), by incorpo- 
rating the conditions Ri = 7Tj (purity of the targets) and 
Ri = R (common mixedness of the sources) in the expres- 
sions for R X R X , T and 5, from Eqs. ©, (|) and (@), 
respectively. These have a particularly appealing form: 



Ry Rx = ^n s n c [n+n_ - 6} , 
t = (i - n c )u + - n s n_ , 



(45a) 
(45b) 



S = y/[(1- K c )Tl + + K s Il-Y - 45K S (1 - TZ C ) , 

(45c) 



where we have defined 1Z C := R 2 cos 2 9, 1Z S := R 2 sin 2 9, 
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FIG. 5: Purifying a pair of mixed states with Bloch vectors of 
length R < 1, separated by an angle 20 G (0, ir] with priorities 
7Ti = 7T2 = 1/2. The plot shows the optimal average fidelity 
for different values of R and the control parameters Si and 
^{1,2,3} f° r R — 0.6. The purification procedure attempts to 
increase, as much as possible, the length R while preserving 
the angle 9. 



(7Ti — 7T2) 2 and 



n 



± 



Tx\ ± 27T17T2 COS 29 . 



(46) 



From the above equations, the indicator function and 
the figure-of-merit can be immediately obtained. At this 
point, though, we specialize to the case S = (i.e., tti = 
7T2 = 1/2) and give explicit formulas in this particular 
case. From Eq. (pi), 



n 



(i - n c )u + - ^n s n c u + u- 

2 [cos 2 9-R 2 cos 9 cos 6 cos (0-6)] , 



(47) 
(48) 



where, in the second line, we used that n + = cos 2 9 and 
II_ = sin 2 9 when m = ir 2 — 1/2. From Eqs. (|l7|) and 



3\ 



HS 



,2 i R 2 sin 2 8 sin 2 8 if O ^ n 

\ + \Rcos(9-9) if SI < 



(49) 

It is now straightforward to see that, if 6 = 6, then 
> with saturation if and only if 6* = 7t/2. That 
is, unless we are trying to purify from antipodal mixed 
states to orthogonal states, the best strategy is always 
a non- unitary transformation (procedure A). The opti- 
mal average fidelity of the purification scheme can be ob- 
tained by using 9 = 9 in Eq. (^). The resulting optimal 
purification performance is shown in Fig. ^| and corre- 
sponds to the best achievable average fidelity allowed by 
quantum mechanics to the purification problem at hand. 

From Fig. M, we see that for small 6, J-hs is typically 
high, regardless of the length R. This can be understood 



in analogy to the fact that collapsing a set of mixed states 
into a single pure state is always perfectly achievable. 
In fact, such a collapse is nearly what is needed in this 
domain, since a pair of pure target states separated by 
a small angle can be well approximated by a single pure 
state. Fig. ^confirms this reasoning by showing that, in 
the small 9 domain, the source Bloch vectors are strongly 
compressed due to the small values of ^{^2,3} and then 
strongly elongated due to the large value of s±. 

For increasing values of 9 the fidelity decreases. Such 
a decay is accentuated if the degree of mixedness of the 
source states is high (small values of R), reflecting the in- 
tuitive idea that it is harder to purify very mixed states. 
In these intermediate regions, a non-trivial combination 
of compressions ^{1.2,3} and translation si of the Bloch 
vectors forms the optimal purifying scheme. Noticeably, 
the optimal procedure has less effect on the qubit (de- 
creasing si and increasing ^{1,2,3}) as 9 increases. 

At 9 = 7r/2, we have = and an optimal uni- 
tary transformation is actually to do nothing (the uni- 
tary transformation V is undone by another unitary V' , 
see Sec. II B 3 ). Note that although the plot of \i\ and /12 
in Fig. g approaches a constant value in between and 
1, the vanishing indicator function introduces a discon- 
tinuity in the purifying operation, since now we should 
use procedure B, hence fix — fi 2 — 1 at 9 = ir/2. Never- 
theless, the values of \x\ and ^2 are utterly irrelevant in 
this case. At this stage both Bloch vectors are aligned 
with the z direction, thus any compression along x and 
y cannot affect the states of interest. 



C. Stabilizing pure states 

A possible use for tracking is to try to cope with the 
presence of noise in quantum computation and commu- 
nication involving qubits. In general, noise processes (we 
restrict ourselves to CP processes) cannot be inverted 
by another CP map, not even when the noise is perfectly 
known p5[ . However, instead of stabilizing the full Bloch 
sphere against noise, one may be interested at stabilizing 
only a limited number of states. Although not perfect, it 
is not uncommon that good stabilization can be achieved 
within this framework. 

In this section, we consider a quantum error correc- 
tion task of this type, which was studied in detail in Ref. 
Jl4| . We will show that the optimal correction scheme is 
merely a particular case of the quantum state purifica- 
tion procedure (with 9^9) introduced in the previous 
section. 

Assume that Alice prepares (with equal probabilities) 
a qubit in one of the non-orthogonal pure states 



l^i) =cos-|+)+sin-|-) , (50a) 
|V> 2 ) = cos^|+)-sm~|-) , (50b) 
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where |±) = ( 1 0) ± 1 1 ) ) / v2 and 9 is the half-angle between 
l^i) and 1^2) in the Bloch sphere representation, hence 
9 £ (0,7r/2). She then sends her qubit to Bob through a 
dephasing channel 



£(0=pS?Z+(l-p)?. 



(51) 



where p is a constant in the range (0, 1/2] that has been 
previously determined. Bob, who does not know which 
of the two states was prepared, has to apply a quantum 
operation so as to ensure that, when Alice performs a 
check-measurement {IV'fcXV'fcl ,1 — iV'fc) (V'fcl} (with k la- 
beling the identity of her actual preparation) on Bob's 
output, the probability of detecting her original prepa- 
ration is as high as possible. This probability equals the 
average fidelity between the possible inputs and the out- 
puts of Bob's operation. 

Our tracking strategy can be of assistance to Bob if 
he regards the two possible noisy states as the source 
states pi = £ (\ipi) and tracks (with equal priorities 
71-j = 1/2) to the target states p 4 = \ipi)(ipi\. In this case, 
the target states are pure and the source states have the 
same degree of mixedness [this follows easily from the 
application of the dephasing map to the states of Eq. 
(p0|)], which is precisely the scenario we considered in 
the last section for quantum state purification. 

The indicator function £1 can then be obtained from 
Eq. ( |48| ) by using the following identities for the angle 
9 (recall that 9 is the half-angle, in the Bloch sphere, 
between the states output by the dephasing noise), 



. „ sin0 (l-2p)cos0 
smO = — and cos0 = , (52) 



R 



R 



where R is the length of the noisy Bloch vectors. Explic- 
itly, 



Q = 2cos 2 9 [l - R 2 + 2psin 2 



(53) 



It is easy to see that, given the ranges 9 £ (0,7r/2) and 
p £ (0,1/2], we have f2 > 0, which implies that Bob 
should always apply the non-unitary procedure A. The 
optimal performance is then obtained by substituting the 
identities (p2|) in the first line of Eq. (E9|) , which gives 



D. Perfectly tracking quantum states 

In this section we evaluate the performance of our 
strategy in circumstances where tracking is known to 
be perfectly achievable. It will be convenient to split 
our analysis in two, namely, the case of two pure target 
states and the remaining cases (in which at least one of 
the target states is mixed). 



1. Pure target states 

In appendix |a] we prove a corollary of Alberti and 
Uhlmann's theorem stating that a CPTP map A per- 
fectly transforming a pair of quantum states pi {i = 1, 2) 
into a pair of pure states exists if and only if pi are 
also pure and 9 > 9. Sin ce our tracking strategy is opti- 



mal (cf. Theorem 111.1), we can infer from Alberti and 
Uhlmann's theorem that it implements tracking with unit 
fidelity whenever R = 1 and 9 > 9. This is explicitly veri- 
fied in the sequence, where the indicator function and the 
figure-of-merit for pure state transformations are com- 
puted. 

We start using Eq. (filf ) with R = 1 (pure source 
condition) to construct the indicator function Q from Eq. 
(^) . After some straightforward manipulation, we obtain 

37Ti7T2 sin cos sin (9 — 9) . (55) 



For our purposes, the only meaningful feature of Q is 
whether it is strictly positive or not, in which case the 
above expression is equivalent to 



n = 2(0-. 



(56) 



since 9,9 £ (0, tt/2] and 7Ti, 1^2 £ (0, 1). Recall that f2 is 
the indicator function obtained in Sec. IIB1, Eq. (Q), 
via an heuristic argument. 

The figure-of-merit, in turn, can be obtained from Eqs. 
(U) and d|) to be = 1 (if Q > 0) and 

1 1 



HS 



27Ti7T 2 cos (29 - 26) (57) 



(if < 0). Note, however, that = 1 if 9 = 6 (i.e., 
Q = 0), in such a way that we can write 



^ HS = ^ + ^/ COs2 ' 



1 - (1 - 2pf cos 2 



(54) 



As expected, this is precisely the optimal fidelity found 
for this problem in . 

It should be clear that our tracking strategy can be 
similarly applied to the stabilization of quantum states 
different from those of Eq. (p0|), prepared with non- 
uniform prior probabilities and undergoing noise dynam- 
ics different from dephasing, in any case still providing 
optimal stabilization. It thus represents a significant ex- 
tension of the results in iH . 



1 if 9 > 9 
if 9 < 9 



(58) 



The first line of Eq. (|58|) is exactly the content of Al- 
berti and Uhlmann's theorem applied to pure state trans- 
formations, whereas the second line establishes the opti- 
mal achievable average fidelity when perfect pure state 
transformation is impossible. 

In conclusion, besides representing a construction of 
Alberti and Uhlmann's map A for perfect pure state 
transformations, our tracking strategy also gives the uni- 
tary map (procedure B) that optimally approximates im- 
possible pure state transformations. 
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FIG. 6: (Color online) Although perfect tracking pi — > 
is physically allowed, the resulting average Hilbert-Schmidt 
inner product [the rhs of Eq. (po|)] is only 0.82 for this trans- 
formation. Our tracking strategy C attempts to maximize this 
number, finding a different transformation which gives an av- 
erage Hilbert-Schmidt inner product [the lhs of Eq. (|6o|)] 
equal to approximately 0.89. As a result, C does not imple- 
ment perfect tracking. 



2. Mixed target states 

The requirement of perfect tracking does not restrict 
the target states to be pure. In fact, the more general 
form of Alberti and Uhlmann's theorem states that for 
any given target states p i; there exists a CPTP map A 
that implements perfect tracking from all source states 
Pi satisfying 



\pi - tp 2 \\tr < \\pi - tp 2 \\tr Vf £ 



(59) 



In contrast to the previous section though, our tracking 
strategy is generally not a construction of the map A in 
this case. As mentioned before, this is a consequence of 
the fact that our figure-of-merit is not as well motivated 
in the case of mixed target states. For example, in sit- 
uations where perfect tracking is possible, the resulting 
average Hilbert-Schmidt inner product does not achieve 
its maximal value. This is further explored next. 

Any CPTP map C implementing perfect tracking must 
satisfy 



7riTr[C(pi)p 1 ]+7r 2 Tr[C(p 2 )P2] = ^i^Pv 



■7r 2 Trp2 



(60) 



Our strategy, though, does not arise from an attempt to 
enforce Eq. (|6Cj|), but instead to maximize its lhs (cf. sec- 
tion II A). Although these actions are equivalent in the 
case of pure target states [the rhs of Eq. ( |60| ) equals 1, 
which is precisely the maximum value of its lhs for states 
satisfying the criterion of Eq. (|S9|)], for mixed target 
states this equivalence is lost. In this case, the lhs can 
typically be made greater than the rhs by employing an 
operation C that elongates the source Bloch vectors to 
nearly pure states, as illustrated in Fig. ||. As a conse- 
quence, the maximization of our figure-of-merit leads to 
a departure from the perfect tracking operation. 

Yet, recall that the average Hilbert-Schmidt inner 
product lower bounds the average fidelity and as such, 



its maximization has some beneficial impact in imple- 
menting tracking, in the sense that it ensures that the 
resulting average fidelity is no less than the maximal av- 
erage Hilbert-Schmidt inner product. 



E. State-dependent Cloner 

One of the most celebrated results in quantum infor- 
mation science is the "no-cloning theorem" ^ || , which 
establishes the impossibility of copying an unknown pure 
quantum state. Since its inception in the literature, a lot 
of work has been done on the topic, both extending its 
range of applicability as well as attempting to weaken 
its impact in practical applications (see [pl| for a re- 
view). Remarkable results in each of these directions are 
the "no-broadcasting theorem" for noncommuting mixed 
quantum states and the Buzek-Hillery optimal quan- 
tum cloning machine 

In this section we consider a state-dependent cloning 
task introduced in Ref. |33| . We will show that our 
tracking strategy provides a straightforward derivation 
of the optimal cloning fidelity obtained in that paper. 
Following §3|, let 



\a) = cos0 |0) + sin0 |1) 
1 6) = sin0|O) +cos0|l) 



(61a) 
(61b) 



for (f> £ [0, 7r/4), be the only two possible preparations of 
a single-qubit, each of which occurring with probability 
1/2. The cloning task is to output the two-qubit state 
\aa) = \a) <S) \a) if the initial preparation is la) or \bb) = 
1 6) (g> 1 6) if the initial preparation is |6). In [ p3[ , a unitary 
transformation U was obtained such that the figure-of- 
merit (the so-called "global fidelity" ) 



1 



(|(aa|C/|aO)| 2 + |(M>|i7|&0)| 5 



(62) 



is maximal. 

The key point that allows the application of our track- 
ing strategy here is that, although the unitary evolution 
U acts on the Hilbert space of a two-qubit system, it was 
shown in B3, Appendix B] that the maximizing U is such 
that U |a0)and U \b0) lie in the two-dimensional subspace 
spanned by {\aa) , \bb)}. Therefore, we can regard this 
cloning as a transformation from the two-dimensional 
subspace spanned by {|a0) , |60)} to the two-dimensional 
subspace spanned by {\aa} , \bb)}. By this same argu- 
ment, we could have even relaxed the condition that the 
system to be cloned is a qubit. 

Let |si) and |s 2 ) ( | tx } and |i 2 )) be the fictitious qubit 
source (target) states, and let 29 (29) be the Bloch sphere 
angle between them. Then, we must have 

(si|s 2 ) = (a0|60) = sin (20) =cos(9, (63a) 
(t x |£ 2 ) = (aa\bb) = sin 2 (2<f>) = cosfl. (63b) 

From the above equations, the angles 9 and 9 can be 
computed in terms of 4>, and the optimal value of F g is 
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given by the optimal fidelity for trackin g between pure 
qubit states, as described in Sec IV D 1 . In particular, 
note that for the present problem, a valid indicator func- 
tion is the one proposed in Eq. (pq), 

= 2 arccos [sin (2(j))} - 2 arccos [sin 2 (20)] < , (64) 

where the inequality holds for the specified range of 
implying that the optimal fidelity is given by Eq. (|5 
with the proper values of 9 and 8, explicitly 



1 



27T17T2 COS f2 . 



(65) 



For 7Ti = 7T2 = 1/2, the above formula can be shown to 
be precisely the same as Eq. (38) of |33| , which gives the 
optimal global fidelity of the doner. Thus we have not 
only reproduced that previous result, but also determined 
how it is optimally modified to incorporate an unequal 
probability of preparation of \a) and \b). 

Finally, let us just mention that the resulting optimal 
tracking unitary operation (call it W) is not quite the op- 
timal cloning unitary operation U appearing in Eq. ( |62| ) 
and detailed in |3j| (U and W do not even act in Hilbert 
spaces of equal dimensions). Instead, W constrains how 
U acts on the states of the form \ipO), but to fully specify 
U we would need to choose U |11) and U 1 01} such that U 
is a unitary matrix. Since this choice is not unique and 
does not affect the fidelity, we can say that W contains 
all the essential information associated with the optimal 
cloning map. 



TRACKING WITH A CONTROL LOOP 



Although the strategy introduced in Section II B has 
been tailored to correspond to a CPTP map, so far no 
insight on how such a map can be physically implemented 
has been given. In this section we provide a realization in 
terms of a quantum control scheme. Namely, procedures 
A and B are shown to have the structure of closed and 
open loop control, respectively. 

We start by giving a possible Kraus decompositions for 
the CPTP maps representing our strategy. This is rele- 
vant here because the Kraus form of a CPTP map enables 
us to interpret that map as some generalized quantum 
measurement (with no record of the outcomes) . For 
n > 0, the transformation C{p t ) = UV(Vp l V jf )LP horn 
procedure A can be written as 



C( Pi ) = (UM^piiUM&y 



(UYM 2 V) Pi (UYM 2 V)i , 
(66) 



with 



Ms . [ ±-1) |+)< + | + Bin (x+!h |_)(_| , 



u-2 = sh, 1 2L_2)|+) <+ |_ 00S (x+2 ]|-)(- 



(67a) 



(67b) 



w here x an d V are defined such th at sinx = M3i cos% = 
yl — /if, sin r\ = p 2 and cos r\ = y/l — p\- 

For fi < 0, the transformation C{pi) = UVpiV^W 
from procedure B is automatically in Kraus form, with a 
single Kraus operator UV. 

We interpret these results as follows. First for fl > 0, 
the unitary V is applied to the system and then a gen- 
eralized quantum measurement with operators M% and 
M 2 is performed. Conditioned on observing the outcome 
'2', a Pauli Y is applied to the system, followed by the 
unitary U. If the outcome is '1', the unitary U is ap- 
plied straight away. Due to this measurement-dependent 
dynamics (feedback), procedure A can be regarded as a 
closed loop control scheme. 

Note that the measurement operators Mi and M 2 are 
not projections, so the implementation of such a mea- 
surement requires the enlargement of the Hilbert space 
(by interaction with an ancilla), with subsequent (projec- 
tive) measurement of the ancilla. Fig. shows a possible 
circuit model for procedure A. 




FIG. 7: A circuit model illustrating the feedback structure 
of procedure A. In the figure, |±i) = (|0) ± i \l))/y/2 are the 
eigenvectors of the Pauli matrix Y, H is the Hadamard gate 
and Zg — exp (—i0Z/2). The highlighted circuit entangles the 
main system with the ancilla and projectively measures the 
ancilla in the basis { , | — i)}. This induces a non-projective 
dynamics of the main system, and for this reason this block 
is referred to as a "weak measurement" . If the measurement 
outcome is then the unitary transformations Y and U 

are applied to the main system; otherwise, only U is applied. 

For £1 < 0, there is clearly no measurement involved, 
hence the control strategy is implemented independent 
of acquiring extra information from the system. For this 
reason, procedure B can be regarded as an open loop 
control scheme. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we have introduced a simple quantum 
version of a common classical control problem named 
tracking. Our quantum tracking problem consists of de- 
termining how to optimally enforce a certain dynamics to 
a qubit system, when the initial preparation of the qubit 
is uncertain (as modeled by a pair of states occurring 
with given prior probabilities) and the desired dynam- 
ics depends on the actual preparation. We presented an 
optimal quantum tracking strategy. 
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The tracking problem studied here is sufficiently gen- 
eral to provide an unifying approach to many prob- 
lems in quantum information science as special cases. 
For example, some cases of quantum state discrimina- 
tion, quantum state purification, stabilization of qubits 
against noise and state-dependent quantum cloning were 
explicitly shown to be instances of quantum tracking. As 
such, previously known quantum limits in the realization 
of these tasks were recovered via the application of our 
tracking strategy. Likewise, our tracking strategy can be 
used to obtain new and improved limits in the realization 
of other impossible quantum machines. 

The derivation of our strategy was largely dependent 
on the fact that our figure-of-merit (the averaged Hilbert- 
Schmidt inner product) is linear in C, which, in turn, is 
constrained to be an element of the convex set of CPTP 
maps acting on qubits. This implies that the optimal 
map C belongs to the subset of extreme points, which has 
been fully characterized in pq| . Thanks to a parametriza- 
tion of these extreme points, the resulting optimization 
problem could be handled analytically when a few mild 
assumptions (supported by numerical observation) were 
made about the form of the optimal solution. The op- 
timality was safeguarded a posteriori via an argument 
based on the SDP structure of the tracking problem. 

Analytical solutions for generalizations of the tracking 
problem studied here (e.g., other figures-of-merit and/or 
larger dimensional quantum systems) seem to require a 
modified approach from the one adopted here. For ex- 
ample, had we chosen to proceed with a better motivated 
figure-of-merit for mixed targets, such as the average fi- 
delity, we would still have the guarantee that the optimal 
C is an extreme point, however optimality results about 
a possible guess would be harder to derive, since it is not 
known if/how the resulting optimization problem can be 
cast as a SDP when source and target states are mixed. 
Alternatively, we could have chosen, for example, to min- 
imize the average trace distance, which can be cast as a 
SDP |55], [jfjj. However, the trace distance is convex in 
C, in which case its minimum is not an extreme point. 
Finally, had we kept our linear figure-of-merit but gen- 
eralized from qubits to qudits for d > 2 (or to multi- 
ple qubits), we would face the problem that the extreme 
points of the set of CPTP maps on higher dimensional 
matrix algebras are not well characterized. 

A possibly simpler generalization is to preserve low di- 
mensionality of the quantum system and linearity in the 
figure-of-merit, but allow for a larger number of possible 
sources and targets. In principle, this problem can be 
approached following exactly the same lines as adopted 
here. In fact, it is not difficult to see that a particular 
case of this more general problem can already be consid- 
ered solved given the results of this paper. Consider we 
are given two sets Si and £2, respectively with n\ and 
712 elements (let N = ni + n-z), of qubit density matrices 
Tj (j = 1, . . . , JV), and want to send every element of Si 
to p.; for i = 1,2. In analogy with Eq. (^]), define the 



figure-of-merit 



^hs - £ Qj Tr [CfoO.fr] + ^ Tr PfaM , (68) 

where the positive numbers qj set the priorities of each 
transformation, and J^jLi 1i = 1- Due to the linearity 
of the trace and of quantum operations, Eq. ( |68| ) can be 
rewritten exactly as Eq. (Q) with 7Ti = X)j=i < 7j' ^2 — 



1 ni 1 
Pi = — 2_, 1j t 3 and p 2 = — V qjTj . (69) 
j=l j=m+i 

Note that 7ti,7T2 > 0, ttl + 7T2 = 1 and p\, p 2 are valid 
density matrices. So, for i = 1, 2 and j = 1, . . . , JV, the 
problem of optimally approximating the JV-state trans- 
formation Si — > /?i with priority qj is equivalent to opti- 
mally approximating the 2-state transformation pi — > 'p i 
with priority 7T,. 
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APPENDIX A: PERFECT TRACKING 
CONDITIONS 

A theorem closely related to the aims of this paper 
has been proved by Alberti and Uhlmann |l3| , consisting 
of a mathematical criterion for the existence of physical 
operations perfectly transforming between pairs of qubit 
states. In this appendix we briefly review this theorem 
and prove an important corollary that is use d in a num ber 
of places in this paper (e.g., sections 



[IB1 and IVD). 



Theorem A.l (Alberti and Uhlmann). Let p\, pi, p~i, 
~p 2 be 2 x 2 density matrices. Then there exists a CPTP 
map A such that 

p 1 =A(pi) and p 2 =A{p 2 ), (Al) 

if and only if 

llPi — *P 2 IUr < ||pi — *Ps|Ur forallteM.+ . (A2) 

where || • ||t r denotes the trace norm. As pointed out 
by Chefles, Jozsa and Winter p7| , the condition ( |A^ ) 
is equivalent to the requirement that the target states 
are no more distinguishable than the source states by 



14 



minimum error probability discrimination (Helstrom p]), 
for any prior probabilities. In the particular case where 
Pi and p 2 are pure states, this just means that the Bloch 
angle between 'p 1 and ~p 2 is smaller than the angle between 
pi and p 2 . This is proved in the following. 

Corollary A. 1.1. Let~p 1 and~p 2 be any two pure distinct 
qubit states separated by an angle 9 6 (0, n] in the Bloch 
representation. Let p\ and p 2 be any (mixed or pure) 
qubit states separated by 9 € (0, 7r]. A CPTP map A 
such that 



p l = A(pi) and p 2 = A{p 2 ) 



(A3) 



exists if and only if p\ and p 2 are also pure and 9 < 9. 



Proof. First note that the inequality (A2) can be equiv- 
alently written with both sides squared. Also, since 
Pi ~ £p2 an d Pi — are hermitian matrices, their trace 
norm can be computed as the sum of their eigenvalues. 
In terms of the Bloch parameters, a straightforward com- 
putation gives 



|Pi-^ 2 || 2 = 4(l 



t 2 - it cos e) 



(A4) 



where we have made use of the fact that t G R + , and 



\pi-tp 2 \\ 2 = 

2 [(1 - t f + (R 2 
+ 2|(l-<) 2 -(i?2-| 



-t 2 R 2 2 - 2ii?ii? 2 cos9)] 

t 2 R\ - 2tR 1 R 2 cose)\ , (A5) 



where Ri gives the magnitude of the Bloch vector for pi, 
i = 1,2. 

Now assume that the absolute value on the right hand 
s ide of Eq. ( |A5| ) can be removed, then the inequality 
(A2) takes the form 



1 + t 2 - 2<cos9 < (1 - t) 2 , 



(A6) 



which for all t E M + is satisfied if and only if cosG = 
1. However, as the (pure) target states are required to 
be distinct, we must have cosO < 1. As a result, the 



inequality (A6) is never satisfied. 

Assume then the complementary case (when the abso- 
lute value of Eq . (|A5 | ) is removed at the cost of a change 
of sign). Then QA2| ) can be written as F(t) < with 

F(t) = (1 - R\)t 2 - 2i(cos9 - R X R 2 cos 9) + (1 - Rf) . 

(A7) 

If R 2 7^ 1, F(t) is a strictly convex function oft, therefore 
cannot be bounded from above by for all t £ R + , so it 
is necessary that R 2 — 1 (p 2 must be pure) . Then, define 
G(t) - F(t)\ R2=11 explicitly 

G{t) = -2i(cos 9 - Ri cos 6) + (1 — R\) , (A8) 

and require G(t) < 0. 

If R\ 7^ 1, G(t) is a linear function of t with strictly 
positive linear coefficient. Again, such a function cannot 
be bounded from above by for all t £ R + , so it is 




FIG. 8: Schematic for proof that i?_ > R y 



necessary to make Ri = 1 (pi must be pure). Finally, 
define H(t) = G(t)\ Ri=1 , i.e., 



H(t) = -2£(cos9 - cos 9) . 



(A9) 



and require H (t) < 0. Clearly, this inequality is satisfied 
for all t £ M + if and only if cos 9 > cos 9, or equivalently, 
9<9. □ 



APPENDIX B: TECHNICAL DETAILS 

1. Properties of 5* and T 

Here, we prove that S > and S + T > if and only 
if one of the following holds 

i) {R±,R 2 } is linearly independent; or 

ii) R 2 } is linearly dependent with T > 0. 
Moreover, we show that the complementary case 

iii) {Hi, R 2 } is linearly dependent with T < 0, 

occurs only if O = 0. 

This result is useful to demonstrate that the coeffi- 
cients /ii, /12, P-3 and si defined in Eq. ( pd] ) for £1 > 
(procedure A) are always (a) well-defined, (b) real (c) 
within the range [0, 1]. We start with the following lemma 

Lemma B.l. Let R\, R 2 be real three dimensional vec- 
tors such that Ri < 1 (z = 1,2). Define R := i?i — R 2 . 
If R- 7^ (i.e., R\, R 2 are distinct), then R 2 _ > R 2 X . 

Proof. Consider the triangle defined by the vectors Ry, 
R 2 and R as shown in Fig. ||. The magnitude of R x 
gives twice the area of the triangle so that 



Ri = h £ R 



2 r>2 



(Bl) 



where h is the altitude relative to the side of length i?_ 
We write the following 



R 2 X = h 2 R 2 _ < mm(Ri,R 2 ) 2 R 2 _ < R 2 _ 



(B2) 



The first inequality is a direct consequence of the 
Pythagorean theorem, and the second follows from R t < 
1 . This establishes that Ri > R\ . This inequality 
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is trivially saturated if i?_ =0. To see that this is 
the only case where saturation occurs, assume i?_ ^ 
and require saturation of both inequalities in Eq. (B2). 



The first inequality is saturated iff R_ — \R{ — R 2 \ 
(by the Pythagorean theorem), and the second one iff 
Rl = R2 = 1. Taken together, these conditions imply 
R- = 0, which contradicts the hypothesis. Therefore, if 
R- ^ (i.e., R x ^ R 2 ), then R 2 _ > R 2 X . □ 

Now, recall that 



S = \J T 2 + AR X (R 2 _ - Rl 



(B3) 



to show that S + R + 7^ 0. Suppose, on the contrary, 

that S = — R,, which is possible only if S = R+ = 0. 
From Eq. (ffl), this can be seen to be equivalent to 
T = R x — R + = 0. To see that this leads to a con- 
tradiction, use once again the fact that R x = R + = 
implies opposing target Bloch vectors of same magnitude 

R, which gives T = —R 2 _R" 7^ 0. The inequality follows 
from the conditions of the problem: the source states 
cannot be identical 7^ 0), and the case where the 
two targets are identical to the maximally mixed states 
has been excluded from the analysis [R 0). 



Assume first linea r ind ependence of {-Ri, -R2} (he., R x 
0). From Lemma B.l, it is immediate that S > 0. More- 



|r|>o, 

(B4) 



S + T = T+yT 2 + AR X (R 2 _ - R x ) > T 



where the first inequality follows from Lemma |B.1| and 
the second is trivial. Therefore, S > and S + T > if 
condition (i) holds. 

For linearly dependent {Ri,R 2 }, it is easy to see that 
S = \T\ and S + T = \T\ + T, therefore S > and 
S + T > if condition (ii) holds. 

To prove the only if part, consider the complementary 
case (iii). It is immediate that S + T = if {Ri, R2} are 
linearly dependent and T < 0, hence (i) and (ii) are the 
only situations where the premise holds. 

It follows trivially from the discussion above that ft = 
for condition (iii) . Simply note that S + T = and the 
linear dependence of the targets Bloch vectors requires 
2R X R X = 0. 



3. Characteristic Polynomials for F 

In this appendix we compute the characteristic poly- 
nomials of the matrix F, Eq. ([33]), with the set of coef- 
ficients given in Eq. (^5|) (for Q > 0, procedure A) and 
Eq. ( p6| ) (for f2 < 0, procedure B). By studying these 
polynomials, we show that F > 0, thus completing the 
proof of the optimality of our tracking strategy. 



a. Procedure A 

For the set of coefficients ( |35| ) (case > 0), the char- 
acteristic equation for F factors as X 2 P2(X) — 0, where 

P 2 (A) = A 2 — r Q A + v [(R 2 _ - R 2 x )Tl - E 2 ] , (B5) 



and 



AR 2 _R X + (S + T) 2 
8RlT 2 a S(S + T) 



(B6) 



2. Well-deflnedness of the dual feasible point 

The propose d va lu es fo r the coeffi cients pi and £3 de- 
fined in Eqs. (35c), (36b) and (36c) have the quantities 
r a and Tf, appearing in the denominator. In this ap- 
pendix we show that this does not lead to any singularity 
as long as the indicated range of fl is considered. 

To see that, note that r a = if and only if R + = 
R x = 0. This, in turn, is equivalent to the statement 
that the targets have Bloch vectors of same magnitude 
R pointing to opposite directions, which used in Eq. (0) 

gives T = —R 2 _R . In these circumstances, ft can be 
easily computed to be = T + \T\ = 0. Therefore, no 



Since both T a and v are positive, the eigenvalues of F 
are n on-n egative if the term in square brackets in the 
Eq. ( p35| ) is non-negative when £1 > 0. We now show 
that this term is non-negative irrespective of the sign of 

n. 

First use Eq. ( J3S| ) to substitute for S 2 , after some 
manipulation we find that 



(R^-R 2 x )Ti-E 2 



R_ 



—2 



/_2 R 2 R 

a R 



S + T 



R 2 _Rl + R 2 ,T 



, (B7) 



where we have defined 



a := 



4i? x 
S + T 



(Ri-Ri 



S -T. 



(B8) 



singularity can occur in Eq. (35c) in the range Q > 0. 

Similarly, r 6 = if and only if T? + = T - 2R X R X , in 

2 

which case we can write fl = S + R + . In the sequence 
— 2 

we show that S + R + > 0, thus no singularity can occur 

in Eqs. ( |36bj ) and (B6c) in the range fl < 0. 

From the definition of S in Eq. (0), it is immediate 
2 

that the inequality S + R, > holds, so we just need first and second summands are non-negative, the term 



with the second equality following from Eq. ua). Note 
that the non-negativity of (i? 2 — R 2 X )T* — E 2 cannot be 
immediately concluded from Eq. (B7) — although the 
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R + T does admit negative values. However, using a 
S — T in Eq. (B7), after some rearrangement we get, 



Rz 



-:2 2 —2 / , . S — T 



SR + + R±R X 1 



S + T 



> , (B9) 



from which the fulfillment of the inequality is obvious. In 
conclusion, procedure A is optimal. 



b. Procedure B 

For the set of coefficients (|36| ) (case Q < 0), the char- 
acteristic equation for F is AP 3 (A) = 0, where 



f 3 (A) = A 3 



and 



(l + RtTi (R 2 _ + Rl)(e 



(RxTt-0 [E*_rlS-R x (?+5*j\ 



8R 4 _r 3 b 



(BIO) 

s 2 ) 

(BH) 
(B12) 



from which it follows that the eigenvalues of F are non 
negative if vd > and u> < tu/ien SI < 0. Next, we 
simplify Eqs. ( |bTT| ) and in order to make it clear 

that these conditions are satisfied. 

It is ju st a matter of applying Eqs. ( j38|) and (|||) to 
Eq. (|Bll|) to show that 



w = i (-n + S + R y R x ) > 0, 



(B13) 



= from which the inequality is clearly seen to hold if SI < 0. 

To prove that u) < if SI < 0, consider first the term in 
the square brackets in Eq. ( B12 ). Again, employing Eqs. 
( p8| ) and ( p9| ) this can be simplified to R^R X T 2 which 
is obviously non-negative. Therefore, the validity of the 
inequality ui < if SI < is now solely conditioned on 
the validity of the inequality 



RxTi~£,>0 for n<o. 



(B14) 



To see that this is so, first note that the only way to 
satisfy the conditions R x = and SI < is to have 
S = T = fl = 0, which implies that (B14) is satisfied 
with saturation. Consider then the complementary case 
fix /0 and SI < 0. Using Eq. (|38| ) for £, and multiplying 
and dividing by 4i? x , we get 



4R> 



—2 



-4i? x i? x (Sl - S) - 4R x Rt 



(B15) 



Now, from Eq. (0), we know that 4 J R x i? 2 _ = S 2 



2 

T 2 + AR 2 X R X , which used in Eq. (B15) gives, after some 
algebra, 



R x r 2 -Z=--^-(S-T + 2R x R x ) >0. (B16) 
. x 

Once again, the inequality is obviously true if SI < 0, 
thus establishing the optimality of procedure B. 
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This is in contrast with the traditional approach to quan- 
tum error correction, which followed the direction of 
adapting classical coding techniques to the quantum do- 
main [ |38j . 

At least not as far as the optimization of the figure-of- 
merit of Eq. (Q) is concerned. A possibly interesting prob- 
lem would be the determination of a figure-of-merit that 
would preserve such behavior for general state transfor- 
mations. 

To see this, note that if the targets are parallel it is 
immediate that SI = \T\ + T and our claims holds 
if T > 0. A straightforward computation shows that 

T = (R 1 + R 2 f- {RTr\ + R%rI + 2R 1 R 2 RiR,2 cos 9), 
which clearly achieves the minimum value iff i2i = 
i?2 = cosB = 1. However, this requires the sources to 
be identi cal, w hich is excluded from the statement of the 



[44] 



[45] 



problem 



[I A 



Therefore, T > 0. 



This follows easily from the inequality R x < R~ (Ap- 
pendix |B~i]) and from R+ = |Si -R 2 \ < - T (the 
equality follows from the anti-parallelism of the target 
Bloch vectors and the inequality from the fact that T < 
for R x = and Q. < 0) 

This follows from the semi-group structure of CP maps. 
An obvious exception arises by restricting to the group of 
unitary noises. In fact, a theorem by Wigner states that 
this is the only exception (see |^] for a proof, see also 

H) 



